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Abstract 

Geometrical optics limit of the Maxwell equations for nonlinear me- 
dia with the Cole-Cole dependence of dielectric function and magnetic 
permeability on the frequency is considered. It is shown that for media 
with slow variation along one axis such a limit gives rise to the disper- 
sionless Veselov-Novikov equation for the refractive index. It is demon- 
strated that the Veselov-Novikov hierarchy is amenable to the quasi- 
classical 5-dressing method. Under more specific requirements for 
the media, one gets the dispersionless Kadomtsev-Petviashvili equa- 
tion. Geometrical optics interpretation of some solutions of the above 
equations is discussed. 

PACS numbers: 02.30.Ik, 42.15.Dp 
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1 Introduction. 

Nonlinear optics provides us with various, very interesting nonlinear phenom- 
ena (see e.g. mi21)- At the same time, it is the source of several nonlinear 
differential equations, modelling such phenomena, with certain very special 
properties. The nonlinear Schrodinger (NLS) equation, which describes self- 
modulation and self-focusing of electromagnetic waves in nonlinear media 
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pi 12] is, perhaps, the most known one. The NLS equation is integrable by 
the inverse scattering transform method and has a number of remarkable 
properties (muhisohton solutions, infinite set of conserved quantities and 
so on). Nowadays there are a number of the, so-called, integrable nonlin- 
ear differential equations which have similar properties [Sill]- A subclass of 
such equations, referred usually as dispersionless integrable equations, has 
attracted a particular interest in the last years (see e.g. |^-[TTj). 

In the present paper we will study integrable structures which arise within 
the geometrical optics of nonlinear media. We will consider nonlinear media 
which are characterized by the Cole- Cole dependence of dielectric function 
and magnetic permeability on the frequency and slow variation of all quan- 
tities along one axis. We will show that the propagation of monochromatic 
electromagnetic waves of high frequency uj in such media is governed by 
the dispersionless Veselov-Novikov (dVN) equation. Namely, we will demon- 
strate that the Maxwell equation in such a situation in the limit uj ^ oo gives 
rise to the standard plane eikonal equation accompanied by the equation 
which describes its deformation in the orthogonal direction. The compati- 
bility of these two equations is equivalent to dVN equation for the refractive 
index. This dVN deformations of the plane eikonal equation preserve, in par- 
ticular, the total "plane" squared refractive index J j n^dxdy. We will show 
also that under more specific conditions, the Maxwell equation is reduced to 
the dispersionless Kadomtsev-Petviashvili (dKP) equation. 

In this paper, we will demonstrate that the plane eikonal equation as 
well as the dVN equation and whole dVN hierarchy, both for complex and 
real-valued refractive indices are treatable by the quasiclassical 9-dressing 
method recently developed in 1131 IT^ . The characterization conditions 
for the 9-data, symmetry constraints for the dVN equation and associated 
hodograph type solutions are also studied. Geometrical optics interpretation 
of some solutions of dKP and dVN equations is discussed too. 

The paper is organized as follows. In section 2 we discuss general proper- 
ties of nonlinear media under consideration. In section 3 the dVN equation 
is derived as the geometrical optics limit of the Maxwell equations. The par- 
ticular case of the dKP equation is discussed in section 4. The quasiclassical 
9-dressing method is applied to the plane eikonal equation in the section 5. 
The dVN equation is treated in the next section 6. Characterization condi- 
tions for 9-data are considered in section 7. Symmetry constraints for dVN 
equation and associated hodograph type solutions are discussed in section 8. 
At the last section 9 we present several exact and numerical solutions both 
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for dVN and dKP equations and corresponding wavefronts. 

Some results of this paper have been announced in the letter [T^ . 



2 Maxwell equations in nonlinear media. 

We begin with the Maxwell equations in the absence of sources (we put the 
velocity of the light in the vacuum c = 1) 

VAH-D = V-D = 

VAE + B = V-B = (2.1) 

and the material equations of a medium 

D = eE, B = /iH. (2.2) 

Here and below V = {d^, dy, dz), while ■ and A denote the scalar and vector 
product respectively. 

We will study a propagation of electromagnetic waves of the fixed, high 
frequency a;, i.e. we will look for the solutions of the Maxwell equations of 
the form ^B] 

E (x, y, z, t) = Eo (x, y, z) e"^-* = Eo (x, z) 

H (x, y, z, t) = Ho (x, y, z) e"^-* = Ho (x, y, z) e--*+-^(-.J/.-) (2.3) 

where Eq, Ho and the phase S{x,y,z) are certain functions and Eq = 
Er=o^~"Eo,„(x,?/,z), Ho = E^^o^"''Ho^„(x,2/,2;). 

Maxwell equations ()2.1ll2.2j) imply the following stationary second order 
equations for Eq and Hq |16| 

V^Eo + uo^neEo + (V log /x) A (V A Eq) + V (Eq • V log e) = 
V^Ho + LU^fieUo + (V log e) A {V A Ho) + V (Ho ■ V log /i) = 0(2.4) 

We assume that the medium is characterized by the Cole- Cole dependence 
[T7j of the dielectric function e and the magnetic permeability /x, namely, 
that 

and 

= /^o + , , ,^ ,2u ^ < < - (2.6) 
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where Eq, /Xg depend on the coordinates x,y,z, while e and ft, depend both 
on the coordinates and the moduh of the fields and their derivatives. In the 
high frequency limit u oo we have 

e = eo{x,y,z) + ^ 
= f^o{x,y,z) + (2.7) 

where ii = (iro)^^''£: and fl^ = (iro)~^''/i. In some textbooks (see e.g. [IHl 
ITU] ) it is argued that at u; — oo the dielectric function has the following 
asymptotic behaviour 

, , const 
e{iu) 1 + — ^, ^ oo. 2.8) 

Experimental results of the Cole and Cole ^\ have demonstrated that the 
theoretical models behind the formula ()2.8|1 are not always correct. The fact 
that the parameter u in the Cole-Cole dependence ()2.5j) and ()2.(jj) belongs to 
the interval (O, |) is of crucial importance for our study. 

At limit u ^ oo, e and /i may depend on the derivatives Sx, Sy, Sz- Differ- 
ent mechanism may be responsible for such a dependence. One is provided 
by a non-local spatial dependence of D on E. Indeed, let us assume that 
a monochromatic wave of frequency uj of the form (j2.3p propagating in a 
medium have a spatial non-local dependence for the displacement D(x, t) on 
E(x, t) of the standard form 



Do(x,c<;)= / / / ^(x'-x,tu)Eo(x'-x,cj)rfx'. (2.9) 




In the particular case where the distribution 5f(x, cj) is proportional to a 
(5-Dirac function 

^(x,cj) =£(x,cj)5(x'-x) (2.10) 

one obtains the usual material equation ()2.2j) . Here, we consider non-local 
effects involving a linear superposition of derivatives 5-functions in the fol- 
lowing way 

oo 

g{^,uj)= J2 ^7^n^,n(x,c.)5(''"^'")(x'-x). (2.11) 

i,m,n=0 
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The distributions S^'"'"^'^^ are defined standardly as 

/(X')5(^'"^'")(X' - X)rfx' = (_l)^+m+n nx 




for a function /(x) defined in and x = {xi,X2,x^). Witfi tlie use of 
the (j2.1ip the formula (j2.9p for the jth component of the field, for j = 1,2,3 
can be written as 

Doj{:si,uj) = ejEoj{:s.,uj) (2.13) 
where the factors ej are given by 



Lm,n=0 



^ , ,\ - \^ ( 1 \^+m+n ^^,m,n(x, iS) g'+"^+"£'oj (x, iS) 



Assuming that in the limit — oo (7000 does not depend on the frequency, 
while ^ 

5(/,m,n(x,U;) » ^2i.+;+m+n ^^.'"."W' ^ OO, (2.15) 

one gets the formula ()2.7|) from ()2.14p . where e and /i depend on S^-, Sy and 

We note that a number of dielectrics has the Cole-Cole dependence of 
dielectric function e on uj pT| EUl E], while for most of studied media the 
magnetic permeability /x does not depend on u. In this paper, however, for 
the sake of completeness we will consider the general case ()2.7|) . Note also 
that at a; ^ oo certain typical effects on nonlinearity, like a second-harmonic 
generation, become irrelevant. 

In addition to the Cole-Cole property we assume that the medium is 
anisotropic one, that is all quantities (E, H, e, /x) vary slowly along the axis 
z, such that, formally, fz = uj~^ f^, (here and below fz = ^ etc.) and ^ is 
a "slow" variable defined by 2; = uj''^. More precisely, we assume that at 
u ^ 00 

fix, y, u^i) = /o(x, y, + u^-'fiix, y, + ^"'V2(x, y,0 + ---- (2.16) 

We would like to note that the Cole- Cole property (j2.7|) is quite generic 
and has been widely discussed, while the particular type of slow variation 
along the axis z given by ()2.16p . seems, did not attract attention before. 
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3 dVN equation as the geometrical optics limit 
of the Maxwell equations. 

Now we will analyze the squared Maxwell equations ()2.4|) for the media de- 
scribed in section 121 i.e. for solutions of the form (j2.3j) . where at uj —>■ oo 

S{x, y, Lu'^O = S{x, y, + ou-^S.^x, y, + u;-^''S2{x, y,0 + ---, (3.1) 
e{x, y, = £o{x, y, + uj-^e^ix, y, + ^^-^-^^(x, + • • • , (3.2) 

/i(x, y, u''"^) = ij,q{x, y, + ^^~>i(a;, v, + c^"^>2(2;, y,0 + ■■■ (3.3) 

Taking into account these expansions, one gets in the leading a;^ order the 
plane eikonal equation 

Si + Si = (3.4) 

where = Eq^q. One, obviously, obtains this plane eikonal equation if there 
is no dependence on z at all. 

In the next 00"^'" and 00"^'"^^ orders, one gets from (j2.4|) the following 
equations 

2SxSix + 2SySiy = EifiQ + Eofii (3.5) 
Si = - {Si + Sfy + 2SxS2x + 2^,52,) + £0/^2 + ^1/^1 + e2^^o■ (3-6) 

As we noted in the previous section, e and depend, in general, on S^.Sy and 
5*2. Having in mind ()3.1ll3.3p . one concludes that in the order uP''^ , Si, Hi 
might depend on the coordinates and only on 5"^. and Sy. So equation ()3.5p 
defines Si in terms of 5*2; and Sy. The function 5*2 in equation ()3.6|) is an 
undetermined function which in turn might depend on Sx,Sy only. Thus, 
equation (j3.6p assumes the form 

Si = ^{x,y,^;Sx,Sy) (3.7) 

where ip is certain function collecting the contributions of 5*1, 5*2 and Ej, fj,j, 
j = 0, 1,2. Let us note that the dependence of the function ip on S is not 
admissible. Indeed, for the solutions of the form ()2.3|) the time-translations 
symmetry t t + const of the Maxwell equations is equivalent to the phase 
displacement S S + const. In order to preserve this symmetry in equa- 
tion (j3.7|) . ip should depend only on derivatives of the phase function 5*. 
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Introducing the complex variables z = x + iy,z = x — iy, one rewrites 
equations (j3.4|) and (j3.7j) as follows 

S,S, = u{z,z,0 (3.8) 
S^ = (p{z,z, ^; 3^,82) (3.9) 

where u = 4n^. 

The condition of compatibility of equations ()3.4|) and ()3.7|) (or ()3.8|) . ()3.9p ) 
imposes constraints on the possible forms of the function if, namely 

S'fV?^ + S'^V^s + MzV^' + MsV^" = M^, (3.10) 

where 

^' = -w^i^^^'^ 3^,82), (p" = {z,z; 3^,82) . (3.11) 

Here we restrict ourself by functions (f which are polynomial in 5*2, Sz and 
compatible with real-valuedness of S and u. For the simplest choice (p = 
ao{z,z,^), equation ()3.10|) obviously gives oq = const, i.e. = 0, S = 
ao^ + S{z, z). For the linear function ip = aSz + aSz + [3 + ]3, one gets 
a = a{z), P = const, and 

u^ = {au)^ + (au)-. (3.12) 

For the quadratic 

^ = aSl + aSl + (3Sz + ^5.- + 7 + 7, (3. 13) 

equations ()3.10p and ()3.8p imply a = 0, /5 = [3{z) and 7 = const, i.e. one 
gets the previous linear case ()3.12p . 
The cubic 

^ = aSl + aSl + pSl + PSI + 7S, + 7^,- + 5 + 5 (3. 14) 

obeys equations ()3.10|) and ()3.8p if 

a = a{z), = 0, 

7g = —UzU — 3aUz, 

6 = const, (3.15) 

and one has the equation 

u^ = {-fu)^ + (ju)-. (3.16) 



7 



In the particular case a = 1 and, consequently 7j = —Su^, equation (j3.16p is 
nothing but the dispersionless Veselov-Novikov (dVN) equation introduced 
in PEj. 

In a similar manner one can construct higher order nonlinear deforma- 
tions of wavefronts along ^-direction which correspond to higher degree poly- 
nomials if. These higher degree cases apparently become physically relevant 
for the phenomena with large values of and 5*^. Thus, if we formally 
admit all possible degrees of Sz and Sz in the right hand side of equa- 
tion ()3.9|) . then one has an infinite family of nonlinear equations, which 
may govern the ^-variations of the wavefronts and "refractive index" u. 
Since equation (|3.9|) should respects the symmetry S ^ —S of the eikonal 
equation (j3.8p . one readily concludes that only polynomials ip of the form 
Sm=i UmS^"^~^ + J2m=i ^mSg"^'^ , are admissible (the constant terms which 
have appeared in the cases ()3.12p . ()3.16|) discussed above are, in fact, irrele- 
vant). Hence, these polynomial deformations are given by 

n 

Si = Yl i^rnS'z^'' + ^mSl^~') > n = 1, 2, 3, . . . (3.17) 

m=l 

where certain functions on u. 

In the case m„ = 1 one gets the dVN equation mentioned above {n = 2) 
and the so-called dVN hierarchy of nonlinear equations. The dVN equation 
has been introduced in [HI as the dispersionless limit of the VN equa- 
tion, which is the 2 + 1-dimensional integrable generalization of the famous 
Korteweg-de-Vries (KdV) equation (see e.g. P 1^1^). 

The dVN equation (|3.16j) and the whole dVN hierarchy have an infinite set 
of integrals of motion The simplest of them is given hj Ci = J J udxdy. 
Indeed, in virtue of equation (j3.16|) 

Ci^ = y j u^dxdy = J j iiiu) ^ + {^u)g) dxdy = (3.18) 

= / (^judz + c.c.) (3.19) 
JdG V2« / 

where G is a domain in C and dG is the boundary of G. For G = C and 
solutions of dVN equation such that 7^ — at |2;| ^ oo 

Cu = j J u^dxdy = 0. (3.20) 
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Thus the dVN hierarchy represents itself the class of deformations of the 
plane eikonal equation (j3.8j) . which preserve the total "plane" squared re- 
fractive index J J n'^{x,y,^)dxdy. The physical meaning of higher conserved 
quantities is not that clear yet. 

Note that for the finite domain G, Ci is proportional to the Dirichlet 
integral over domain G, Ci = 4 f J (5*^ + Sy) dxdy. The formula (j3.18|) 
gives us the variations of the Dirichlet integral due to the dVN deformations. 

4 Integrable deformations of the quasiplane 
wavefronts via dKP equation. 

Let us consider a more specific situation in which the propagation of electro- 
magnetic waves in the media discussed above exhibits also a slow variation 
along the axis y, namely dy = edrj, where ?7 is a slow variable defined hy y = - 
and e is a small parameter. Let the phase function 5* and u in the eikonal 
equation (j3.4p have the following behaviour as e ^ 

S (x,^,^^ = ^-S (X, 7], , 

{x,^,{j = ^ -q{x,r],^) . (4.1) 

Considering the simplest cubic case ()3.14|1 for a = 1 in Cartesian coordinates 

= ^Sl - ^S^S'y + UA + U,Sy (4.2) 
where 7 = f/i + iU2, and assuming that 

^1 = tIt + 7^' = -2^^' (4-3) 

loe^ 4 

one gets from (j3.4|) and (|4.2p the following equations 

S, = Sl + q, (4.4) 
= Sl + ^qS, + U. (4.5) 

Compatibility of equations ()4.4|) and ()4.5p gives rise to the well known dis- 
persionless Kadomtsev-Petviashvili (dKP) equation = ^qq^. + |5^^gr?r? and 
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U = ^d'^Qy. The dKP equation is rather well studied (see e.g. [H]-[TI1] 
and reference therein). The KP equation itself represents the most known 
2+1-dimensional integrable generalization of the KdV equation. 

In the general case, assuming that the polynomial ip (x, ^, ^] Sx, eSn) has 
an appropriate behavior as e — > 0, one readily shows that equation ()3.7|) is 
reduced to 

S^ = ip(^x,r],^]S^^ (4.6) 

where (p is an odd order polynomial in S^- The condition of compatibility 
between equation (j4.4|) and equations of the type (j4.6p gives rise to entire 
dKP hierarchy. These equations describe propagation of the quasi-plane 
wavefronts y = const in a medium with "very large" refractive index. 

The dVN and the dKP equations being relevant in particular situations 
of propagations of waves in certain nonlinear media have an advantage to be 
integrable. 



5 Quasiclassical t^-dressing method for the plane 
eikonal equation. 

In this and next sections we will demonstrate that the plane eikonal equation 
and dVN hierarchy both for complex and real refractive indices are treatable 
by the quasiclassical (9-dressing method. 

The quasiclassical (9-dressing method is based on the nonlinear Beltrami 
equation [T^-[Ti] 

S-, = W{\,\;Sx) (5.1) 

where ^(A, A) is a complex valued function, A is the complex variable, S\ = 
II and W (the quasiclassical (9-data) is an analytic function of 

oo 

W (A, A, Sx) = J2 "^niK A) {S^r , (5.2) 

n=0 

with some, in general, arbitrary functions w„(A, A). 

To construct integrable equations one has to specify the domain G (in 
the complex plane C) of support for the function W (W = 0, X E C/G) 
and look for solution of (jS.ip in the form S = Sq + S, where the function 
5*0 is analytic inside G, while S is analytic outside G [121- [EI. order 
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to construct the eikonal equation on the plane, we choose G as the ring 
P = {A G C : i < |A| < a}, where a is an arbitrary real number (a > 1), 
and select B — type solutions satisfying the constraint 

S[-\-\) =-S{X,\) . (5.3) 

Then, we choose 

^0 = + ^. (5.4) 
A 

Due to the analyticity of S outside the ring and the property ()5.3p one has 

oo ^{oo) 

S = ^^^\ A->oo, (5.5) 

n=0 

and 

oo 
n=0 

In particular, S (0, 0) = 0. 

An important property of the nonlinear (9-problem 1)5.11) is that the deriva- 
tives f = St of S with respect to any independent variable t, obeys the linear 
Beltrami equation 

h = W'{X,X;Sx)fx (5.7) 
where W (A, A; 0) = ^ (A, A; 0). Equations ()5.7j) has two basic properties, 
namely, 1) any differentiable function of solutions /i, . . . , /„ is again a solu- 
tion; 2) under certain mild conditions on W, a bounded solution / which is 
equal to zero at certain point Aq € C, vanishes identically (Vekua's theorem) 



These two properties allows us to construct an equation of the form 
^(ySzjSz) = 0, with certain function Q. Indeed, taking into account ()5.4|1 . 
one has 

Sz = \ + S,, (5.8) 

Sz = \ + S,, (5.9) 

i.e. Sz has a pole at A = oo, while has a pole at A = 0. The product SzSz 
is again a solution of the linear Beltrami equation ()5.7p and it is bounded on 
the complex plane since 

SzSz = l + jSz + XSz + SzSz (5.10) 
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where Sg {X = 0) = and XSz —>■ S[^^ as X ^ oo. 

Subtracting 1 + s[^^ from the r.h.s of equation ()5.1Up . one gets a solution 
of equation ()5.7|) which is bounded in C and vanishes as A — >• oo. According 
to the Vekua's theorem it is equal to zero for all A. Thus we get the equation 

SzSz = u{z,z) (5.11) 

where 

u = l + S[f. (5.12) 

In the Cartesian coordinates x, y defined hj z = x + iy, equation ()5.1ip is 
the standard two-dimensional eikonal equation 

(V5)' = (5.13) 

where V = |) and = f . 

Using the 5-problem 1)5.11) . one can, in principle, construct solutions of 
equation ()5.1H) . So, the quasiclassical (9-dressing method allows us to treat 
the plane eikonal equation ()5.11)1 in a way similar to dKP and d2DTL equa- 
tions [121- [EI- laote that the phase function S in ()5.1H) depends also on 
the complex variables A and A. Curves S{z, z) = const define wavefronts. 
The 9-dressing approach provides us also with the equation of light rays. 
Indeed, since r.h.s. of fl5.13|) does not depend on A and A, the differentiation 
of ()5.13|) with respect of A (or A) gives 

V5-V0 = O (5.14) 

where (j) = Sx {or (f) = Si)- So, the curves S = const and = const 
are reciprocally orthogonal and, hence, the latter ones are nothing but the 
trajectories of propagating light. Thus, the 9-dressing approach provides us 
with all characteristics of the propagating light on the plane. Note that any 
differentiable function {Sx, Sx) is the solution of equation ()5.14|) too. 

Typically the refractive index is real and positive one being defined in 
terms of high frequency limit of dielectric function and magnetic permeability 
[TBj . However, it was noted in i24j that for certain media the product eo/^o 
and, hence, u = can be negative as well. Moreover, one of the ways to 
describe the damping effects is to consider a complex- valued jTHllSlI- So, 
models of optical phenomena deal both with real- valued and complex-valued 
refractive index. 
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In general, within the 9-dressing approach one has a complex- valued 
phase function S and, consequently, a complex refractive index. To guar- 
antee the reality of u, it is sufficient to impose the following constraint on S 

S{X,X)=s(^j,jy (5.15) 

Indeed, taking the complex conjugation of equation 1)5.111) . using the differ- 
ential consequences (with respect to z and z), of the above constraint and 
taking into account the independence of the l.h.s. of equation ()5.1H) on A, A, 

one gets 

u {xn) = (A, A) 5*^ (A, A) = S2 Q, ^ 5*2 (^j, jj = u (xn) , (5.16) 

i.e. the "refractive index" u is real one. The constraint ()5.15j) leads to the 
relations ^2n+i = Moreover, this constraint implies also that the 

function S is real-valued on the unit circle |A| = 1 (S* (A, A) = ^(A, A), 
|A| = 1). This provides us with the physical wavefronts. 

The 9-approach reveals also the connection between geometrical optics 
and the theory of the, so-called, quasiconformal mappings on the plane. 
Quasi-conformal mappings represents themselves a very natural and impor- 
tant extension of the well-known conformal mappings (see e.g. [2^1 12^1 )• In 
contrast to the conformal mappings the quasi-conformal mappings are given 
by non-analytic functions, in particular, by solutions of the Beltrami equa- 
tion. 

According to [23 1211 ^ solution of the nonlinear Beltrami equation ()5.ip 
defines a quasi-conformal mapping of the complex plane C. In our case we 
have a mapping S (A, A) which is conformal outside the ring G and quasi- 
conformal inside G. Such mapping referred as the conformal mapping with 
quasi-conformal extension. So, the quasi-conformal mappings of this type 
which obey, in addition, the properties ()5.3|) . ()5.4|) and ()5.15p provide us with 
the solutions of the plane eikonal equation ()5.1ip . In particular, wavefronts 
given by ^(A,^;^,^) = const are level sets of such quasi-conformal map- 
pings. In more details, the interconnection between quasi-conformal map- 
pings and geometrical optics on the plane will be discussed elsewhere. 
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6 Quasiclassical c^-dressing method for dVN 
hierarchy. 

In this section we will apply the quasiclassical 9-dressing method to the dVN 
hierarchy. For this purpose we consider again the nonlinear 9-problem ()5.1|) . 
fl5.2p on the ring G with the constraints ()5.3|) and ()5.15|) and introduce inde- 
pendent variables Xn, Xn via 

oo oo 

So = Y, x„A2"-^ + J] x„A-2"+^ (6.1) 

n=l n=l 

The derivatives S^^ and S^^ obey the linear Beltrami equation ()5.7|) . and 
using the Vekua's theorem one can construct an infinite set of equations of 
the form 

Q Xn, Sx„, Sx„) = 0. (6-2) 

Repeating the procedure described in the previous section, one gets the plane 
eikonal equation for z = Xi. In a similar manner, for the variable ^ = X2 = X2, 
taking into account that 5*^ = A^ + ^ + S^, one obtains the equation 

= Si + Si + VS, + VS, (6.3) 

where V = — 35*1^^ = —Sd^^Uz- Evaluating the terms of the order A~^ in 
the both sides of equation ()6.3|1 . one gets the dVN equation 

= -3 {ud^'uz)^ - 3 {ud:'u,).^ . (6.4) 

Considering the higher variables ^„ = x„ = x„, (n = 2, 3, . . . ), one constructs 
all equations ()3.17p . and hence, the whole dVN hierarchy. It is a straight- 
forward check that the constraint ()5.15p is compatible with equations ()6.3p 
and ()3.17p . So, the 9-dressing scheme under this constraint provides us with 
the real- valued solutions of the dVN hierarchy. 

If one relaxes the reality condition for u (i.e. does not impose the con- 
straint (j5.15p ). then there is a wider family of integrable deformations of the 
plane eikonal equation with complex refractive index. To build these defor- 
mations one again considers the 9-problem ()5.1|) . chooses the domain G as 
the ring (as before), imposes the S-constraint ()5.3p . but now chooses 5*0 as 
follows 

oo oo 
n=l n=l 
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where yi = Xi = z. Repeating the previous construction, one gets again the 
eikonal equation (|5.1ip . but now with a complex- valued u. Considering the 
derivatives Sx„, Sy„ with n = 2,3, . . . , one obtains the two set of equations 

n 

Sxr^ = '^Urn{Xn,yn) (S^f"'''^ , n = 2,3,... (6.6) 

m=l 

n 

Syn = ^^miXn^Vn) (S^f""'^ , 71 = 2,3,... (6.7) 



m=l 



Equations (|5.1ip and (j6.6j) give rise to the hierarchy of equations 

Ux„ = Fn {u, U^, Uz) , (6.8) 

the simplest of which is of the form 

u^^ = {udj^u^) ^ (6.9) 

where 

Sx, = Sl + VS,. (6.10) 

Equations (jS.llj) and ()6.7j) generate the hierarchy of equations Uy^ = F„, the 
simplest of which is given by 



u. 



y2 



-3{ud:'u,).^, (6.11) 



and 

Sy, = Sl + VS-,. (6.12) 

For both of these hierarchies the quantity J udz A dz is the integral of 
motion as for the dVN hierarchy. It is easy to see that equations ()6.9|] 
and ()6.11|) imply the dVN equation ()6.4|) for the variable ^ = ^^j^. 



7 Characterization of (9-data. 

As we have seen, the constraints ()5.3p and (jS.lSj) guarantee that one will 
get the eikonal equation ()5.1H) with real valued u. In this section we will 
discuss the characterization conditions for 9-data W (A, A, S\j which provide 
us with such result. In general, if one considers 5-problem with a kernel 
defined in a ring and the function 5*0 singular at two points, e.g. A = 
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and A = oo, without 5-constraints (|5.3p . one constructs the dispersionless 
Laplace hierarchy ^3] associated with the quasiclassical hmit of the Laplace 
equation, i.e. with the equation 

Sz{X, z, z)Sz{X, z, z) — aSz = u{z,z), VA G C (7.1) 

where a = dzS{0,0) and u = 1 + dzS[°°\ The eikonal equation (j3.4|) is 
obtained as a reduction of equation ()7.1|) taking 5'(0, 0) independent on z. 
In fact, the B condition ()5.3|) . producing 5(0,0) = realizes this reduction. 
In what follows we will discuss how one has to choose the d-data. in a way 
to construct the dVN hierarchy directly. We will find constraints which 
are dispersionless analog of the constraints found in ^7, , which specify two- 
dimensional Schrodinger equation with real-valued potential. In particular, 
we will see that it is possible to weaken slightly the B-condition ()5.3|) . since 
the value S{0, 0) is fixed up to a constant by dVN reduction. 

In the following we will focus on solutions of d problem of the form S = 
Sq + S, where Sq has polynomial singularities at A = and A = oo and S is 
holomorphic at these points and such that 

lim 5 (A, A) =0, (7.2) 

lim S (A, A) = 5(0, 0, z, z) = -iv{z, z), (7.3) 

A=0 

Lemma 1 Let the kernel W in equation / I5.i)) satisfies the assumptions of 
the Vekua's theorem, and let S be its solution. The condition 

S {X,X) = -S {-X,-X) + const (7.4) 

is verified if and only if 

W (A, A, Sx (A, X))=W' (-A, -A, -Sx (-A, -A)) (7.5) 



Proof. The condition ()7.5|) is a necessary one. Indeed, starting from fIS.ip . 
one obtains 



^ {Sz (A, A)) = W (A, A, Sx (A, A)) — {Sz (A, A)) (7.6) 
t [Sz (-A, -A)) = W (-A, -A, -S, (-A, -A)) |- [Sz (-A, -A)) (7.7) 
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Exploiting the condition (j7.4j) in (j7.7p . one gets the equahty (j7.5|) . 

The condition (j7.5p is sufficient. Indeed, let us introduce the function 

$ (A, A, z, z) = (A, A) + (-A, -A) . (7.8) 

Using the equations ()7.2|1 and ()7.Hj) . one has 

lim $ (A, A, z, ^) = 0, (7.9) 
\im<l> {X,X,z,z) = -2iv,{z,z) . (7.10) 

A— >0 

Both terms in the right hand side of (j7.8p satisfy the Beltrami equation (j7.6|) 
and ()7.7|) respectively, from which, exploiting the constraint ()7.5|) . one con- 
cludes that 

H (A, A, z, z) = W (A, A, Sx (A, A)) ^ (A, A, z, z) . (7.11) 

The function $ is a solution of the Beltrami equation vanishing at A ^ oo. 
So, $ vanishes identically on whole A-plane, so 

f (U)--f (7.12) 

In particular $ (0, 0, z, z) = —2ivz {z, z) = ,that is Vz {z, z) = 0. Analo- 
gously it is possible to demonstrate that 

H.-A), (7,13) 

Hence Vz {z,z) = 0. Equations ()7.12|) and ()7.13|) lead to the relation ()7.4p □ 

Lemma 2 Let S (A, A, z, z^ be a solution of d problem i5.1}) such that \7.^ 
and ([7.4 ^ are verified. Then the condition 

S{X,X)=s(l,j) +tviz,z) (7.14) 



is verified if and only if 



X'W (A, A, Sx (A, X))=-W (^i, J, -X'S-, Q, ^ j ^ . (7.15) 
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Proof. The condition ()7.15j) is a necessary one. Let us consider the 
complex conjugation of equation ()5.ip 

dx {S (A, \))=W (A, A, S, (A, A)) . (7.16) 

Since 

S (A, A) = 5 + tv{z, z) = S (e, I) + iv{z, z) (7.17) 

where ^ = A ^ and \ = A^^, the left hand side of ()7.16j) can be written as 
follows 



1 /I 1 .2. ^ /I 1 

A^'^ VA' A' "-^^ VA'A 



^W{-,-,-yd-xS(-,-)), (7.18) 



that provides us with equation ()7.14|) . 

The condition ()7.15|1 is a sufficient one. The 5-equation ()5.ip written in 
terms of the variables ^ = A ^ and ^ = A~^ 

S-^ = W{^,lS^{^,l)) (7.19) 

is equivalent to 

Using equation ()7.16|) . multiplied by A^, and equation ()7.15p . one concludes 
that 

dxS {X,X) = dxS . (7.21) 

Integrating ()7.21|) . one gets 

S {xrx) = S +v{z,z). (7.22) 

Let us note that v cannot depend on A since the function S is meromorphic 
outside the ring. Now, evaluating the equality ()7.22|) at A = and A — > cxo, 
one obtains _ 

V {z, z) = S{0, 0, z, z) = -S{0, 0, z, z). (7.23) 
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So, f is a purely imaginary function 

v{z,z) =iv{z,z) (7.24) 

where f is a real valued function. This complete the proof □. 

Note that Vz = ivz{z, z) = —a, in other words, it is the coefficient in front 
of the "magnetic" term in equation ()7.1|) . When it vanishes (i.e. a = 0), one 
has the pure potential equation ()7.1|) . that is the eikonal equation. 

Combining together the lemmas (0) and one gets the following theo- 
rem: 

Theorem // the d-data W of the d-equation \5.1\) obey the constraints 

AW(A,A,5a(A,A)) =-iy Q,^,-A'5xQ,^^^ , (7.25) 

W (A, A, Sx (A, A)) = W (-A, -A, -Sx (-A, -A)) , (7.26) 

then this d-problem provides us with the eikonal equation with real-valued 
refractive index. 



8 Symmetries, symmetry constraints and re- 
duction method for dVN equation. 

The dVN equation and the dVN hierarchy possess a number of remark- 
able properties typical for the dispersionless integrable equations and hier- 
archies ITK] . Here we will discuss some aspects of symmetry proper- 
ties of the complex dVN equation, symmetry constraints and corresponding 
1 + 1-dimensional equations which provides us with solutions of the eikonal 
equation and the dVN equation. 

We will concentrate on equation ()6.9|) . i.e. equation 

= -3 {ud^^u,)^, 

with the complex-valued u, which is the compatibility condition for the sys- 
tem 

SzS-z = u, (8.1) 
S^, = Si + VSz (8.2) 
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where Vg = —3uz- Infinitesimal continuous symmetries u ^ u + Su of equa- 
tion (j6.9|) are defined, as usual, by its linearized version, i.e. by the equation 

C5u = (8.3) 
where the linear operator C acts as follows 

C = 9,, + 3d, {{d^'uz) (■)) + 39, {ud^'d^i-)) . 

The system fl8.HI8.2p is quite relevant for an analysis of equation ()8.3|) . 
Namely, it is straightforward to check that a class of its solutions is given by 

TV 

6u = y^cJs,-S,) (8.4) 

i=l 

where Si, Si are arbitrary solutions of the system ()8.1|) . ()8.2|) with given u, 
Ci are arbitrary constants and is an arbitrary integer. 

The formula (j8.4p provides us with a wide class of symmetries of the 
dVN equation. In particular, one can choose Si = S {X = Aj), and Si = 

S (^\ = Xij. In the case Aj = Aj + /Xj, where /Zj ^ and q = ^, one has a 

class of symmetries given by 

N 

Su = ^Ci(f)i,, (8.5) 

i=l 

where 0j = ff ('^) |;s^_;^.- the simplest case n = 1 one has 

Su = ccj).,^, (8.6) 

i.e. a symmetry is given by the divergence of the vector {(j)x,(l)y), which is 
normal to the light rays (p = const, or tangent to the wavefronts S = const. 

Note that one has the same class of symmetries also for the real dVN 
equation ()6.4p . 

Any linear superposition of infinitesimal symmetries is obviously a sym- 
metry too. In virtue of (j8.3p . the requirement that certain linear superposi- 
tion of symmetries vanishes, is apparently compatible with the dVN equation. 
An obvious set of symmetries of the dVN equation is given by Ux„, Uy„. Com- 
bining them with the symmetries ()8.4p . one gets the following set of possible 
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symmetry constraints 

N 

{anUx„ + PnUyJ = ^ (^Si - ^_ (8.7) 
n 1=1 

where a„ and /3„ are arbitrary constants. Symmetry constraints ()8.7p are 
compatible with the dVN hierarchy and reduce it to the famihes of the 1 + 1- 
dimensional equations similar to the dKP and 2DTL case |28j . 

Here we will consider the simplest symmetry constraint of the form 

= 0,,. (8.8) 

Assuming that the integration "constant" is equal to zero, one has 

u = 0,. (8.9) 

Equations ()8.H) and ()8.2|) imply 

S,S, = 0„ (8.10) 
0,^2 + ^,02 = 0, (8.11) 
3^,V. + ^0. (8.12) 



and 

V = -30,. (8.13) 
Substitution of into ^TH^ gives 

0, = -Si (8.14) 

The compatibility condition for equations (|8.10|) and (|8.14p . readily, leads to 
the equation 

S,Sl = (8.15) 
where a is an arbitrary constant. So 

S, = aS;\ (8.16) 

while equation ()8.10j) becomes 

0_ = m-|0j. (8.17) 
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Equations (|8.2j) and (j8.12j) take the form 



S.,=ASl (8.18) 
K = -60'- (8.19) 

It is easy to check that equations ()8.17|) and ()8.19p are satisfied in virtue of 
equations ()8.16|) . ()8.18p and the relation ()8.14|) . Thus, we have proved that 
any common solution S{z,z,X2) of equations 

S, = aS'J (8.20) 
S., = ^Sl, (8.21) 

provides us with a solution of the complex dVN equation (|6.9|) . given by 

u = aS'i and V = 35*^. Of course, any solution S{z,z) of equation ()8.20p 
gives us a solution of the eikonal equation ()8.H) with the refractive index 

2 

u = aS'i . 

One can rewrite equations ()8.2()j) and ()8.2H) in the form of quasilinear 

2 

equations for ij) = oi directly in terms of m = aSi . One has 

u-z = -^u-\,, (8.22) 
u^^ = 12a-^u\,. (8.23) 

It is an easy check that a function u which obeys both these equations is a 
solution of equation ()fj.9|l . Equations ()8.22j) and ()8.23j) are 1 + 1-dimensional 
quasilinear equations. So, the problem of construction of solutions of the 
2 + 1-dimensional dVN equation is reduced to the construction of common 
solutions of 1 + 1-dimensional equations. Such a method of construction of 
solutions for dispersionless integrable equations is known as the reduction 

method lEiiziiniEHiEniisiiiisiiini. 

The method of characteristics provides us with general solution of equa- 
tions ()8.22|) and ()8.23p in the hodograph form 

Q,3 

z u-^z + 12a-^u^X2-ip{u) = (8.24) 

3 

where ip{u) is an arbitrary function. Solutions of equation ()8.24|) give us a 
class of solutions of the complex dVN equation (|6.9|) . In a similar manner 
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one can treat the symmetry constraint Uz = (p^zi equation (|6.1H) . Namely, 
one gets equations (|8.2(JII8.2^ with the substitution z z and X2 ^ y2- 
Under more general constraint 



N 

^c,(^,-^,)^^_, (8.25) 



one has 



« = (8.26) 

i=l ^ 
N 

V = -3J2c^{s,-Si^^, (8.27) 
and equations (j8.ip and at n = 2 are reduced to the following systems 

N 

SkzSkz = Ci (^Si — Si^ _ , (8.28) 

i=l 
N 

SkzSk-z = c^ - , k = l,...,N, (8.29) 



i=l 



and 

N 

Skx2 = — •iSkz ^ Ci (^Si — Si^ (8.30) 

i=l ^ 
N 

Sux, = Si - 3Skz - ^i) ^ k = l,...,N. (8.31) 

i=l ^ 

Common solutions of equations 1)8.28118.3111 provide us with solutions of the 
complex dVN equation given by ()8.2(ip . 

These type of systems, symmetry constraints for the real dVN equa- 
tion ()6.4p and associated constructions within the quasiclassical (9-approach 
will be discussed in the paper 
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9 Examples of solutions: wave fronts and re- 
fractive indices. 



In this section we shall present, as illustrative examples, several exact and 
numerical solutions of the eikonal equation, the dVN and dKP equations 
and visualize wavefronts distribution and refractive indices. For this end, we 
will use the solutions of dKP hierarchy obtained by quasiclassical 9-dressing 
approach in jMj, and solutions of dVN hierarchy obtained in the previous 
section. 



9.1 dKP wavefronts. 

In order to interpret the S function associated with the dKP equation as a 
perturbation of a plane wavefront, it is convenient to rewrite equations ()4.ip 
as follows 

f3S = r]-f3S, (9.1) 
Pn' = ^-Pq (9.2) 

where f3 = 2e^ is a small parameter, and S satisfies equations ()4.4j) and ()4.5|) . 
In what follows we will discuss a two-dimensional explicit solution and its 
three-dimensional generalization. In our interpretation, the two-dimensional 
solution describes wavefronts which are plane, that is, they do not depend on 
^ variable. The three-dimensional solutions show us slow deformations along 
^-axis. 

The purely two-dimensional case is given by 

S = -^-^^-c; |A|<1, (9.3) 
2 a - 2A II-' K J 

g = - (4x + 16?7^) (9.4) 



1 , X (x + 4?7^) 

a=-, b = 2v--, c= ^ " . 9.5 
7] zr] At] 



where 



It is associated with the solution obtained in the paper [Hlj for the particular 
choice of the kernel W = Sj^. Calculating S" at A = 0, we obtain a real and 
non-singular phase function on the x, ?7-plane 

^(A = 0) = -AxT]. (9.6) 
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Figure 9.1: Wavefronts distribution for different values of ^ = 

0.0,0.04,0.08,-0.04,-0.08. For a clearer visualization the variable r = lO^ry has 
been used. The shady region corresponds to the strip < for the figure a) and 
\P{S.)x + Q{£,W\ < ^ for the figure b). 



The equation of the corresponding wavefronts is the following 



X 



k — 7] 
4(37] 



(9.7) 



where k is an arbitrary constant, whose range of validity is restricted to the 
region where ^ ^ 1. In this case it is the strip |x| <^ The figure IHUl-a) 

shows the distribution of wavefronts (obtained setting different values of the 
parameter k), and the shady region represent the strip |x| < where the 
approximation reasonably holds. The figure 12121 -b) represents the refractive 
index density. Light rays turn off towards darker regions corresponding to 
larger refractive index. 

The three-dimensional generalization of the previous case has been dis- 
cussed in the same paper jSlI. In particular, from the formula (35) one gets 
the following real phase function at A = 



A=0 



27] {PiOx + QiOv') 



(9.8) 
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where 



3^0 V(go-4po)^ |go|(go-4po) 
p, = i + 36e-(l-120^ go = 3(l + 40. 

Note that for = the expression for S ()9.8|) coincides with expression ()9.6|) . 
Moreover,the corresponding dKP solution is 

4 fs - 12^ + 4 (1 - 12^5^ + X - Ar]^) 

n - 1 ^ . (9.9) 



3(1+40 (l - 12^ + 2(1 - 120 



Of course, expression ()9.9j) becomes ()9.4j) when ,^ = 0. We discuss the ^- 
deformations of wavefronts in a sufficiently small neighbourhood of the plane 
^ = 0, in order to avoid ^-values for which q blows up or becomes complex 
valued. 

The explicit expression of the wavefronts is the following 

_ , + 2miW - k 

where k is an arbitrary constant. The wavefronts ()9.10|) should be considered 
inside a strip given by 

\P{0^ + QiOv'\<^^- (9.11) 

The figure I^D shows the wavefronts distribution inside a reasonable strip 
for different ^-values. The figure l^^ demonstrates that the variation of refrac- 
tive index is more significant for ^ > 0, where we observe a larger deformation 
of the wavefronts. One can see this in the figures 19.31 and 19.41 



9.2 dVN wavefronts. 

Here, we analyze a refractive index obtained by solution of the hodograph 
relation (j8.24|) . and a corresponding numerical solution of eikonal equation. 
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Figure 9.2: Density plot comparison of the refractive index u for ^ = —0.08,0.0,0.08, 
where r = lO^rj. Darker regions correspond to a larger refractive index. 



-0.08 




Figure 9.3: Deformation of a wavefront at different values of ^. 




Figure 9.4: Deformation of a wavefront along the t = direction, where r = lO^t]. The 
wavefront bends in such a way that light rays turn off towards regions where the refractive 
index increases. 
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-10 -5 5 10 -10 -5 5 10 

X X 

Figure 9.5: Density plot of the a) real and b) imaginary parts of the refractive index n 
on the plane X2 = 0. 

Note that equation ()8.24|) provides us with the complex refractive index and, 
by consequence, the complex phase functions describing, in a way, damp- 
ing effects of the electromagnetic wave (see equation ()2.H|l ). Assuming m 7^ 
and choosing, for simplicity, = 3, we have the following form of equa- 
tion (jH:^ 

-X + iy + [x + iy)u^ + 4x2^^ - u^ip{u) = 0. (9.12) 

A particular choice of the function ■?/'(«) sets the behaviour of u on the bound- 
ary X2 = 0. Assuming ip polynomial on u, the numerical analysis shows qual- 
itatively similar behaviours for different polynomial degrees. As particular 
example, we set = 2m + ^u^. The figure 1^31 shows the density plot of the 
real and imaginary parts of a refractive index n = | on the transverse section 
X2 = , obtained by numerical solution of ()9.12p . Representing the phase 
function in terms of its real and imaginary parts, S = Si + iS'2, one rewrites 
the eikonal equation (|8.1|) as the system 

SlxS2x + SlyS2y = 2U2 (9.13) 

where u = Ui +iu2- A numerical solution of the system (j9.13p for the refrac- 
tive index plotted in the figure 19.51 It has been obtained in the range x G 
[—10, 10] and y G [—10, 10] exploiting the function NDSolve of the Standard 
Mathematica Packages with boundarv conditions 5*1 flO. v) = 82(10, y) = 
y^ and Si{x, —10) = Si{x, 10), S2{x, —10) = S2{x, 10). The figures IHIHl-a) and 
19.61 -b) shows 5*1 and S2 respectively. In particular, one should to note that the 
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10 



10 



Figure 9.6: a) Real part of the phase function with boundary conditions 5i(10, y) = 
and <S'i(a;, —10) = Si{x, 10). b) Imaginary part of phase function with boundary conditions 
52(10,2/) = y2 and S2{x, -10) = S2{x, 10). 




regions where 5*2 increases, correspond to the damping of the electromagnetic 
wave. 

Finally the figures 19.71 -aVb) and c) show the wavefronts configuration at 
different values X2 = 0, 5, 10. 
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